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Abstract In the first sections of this article, we discuss two vari- 
ations on Maxwell's equations that have been introduced in ear- 
lier work — a class of nonlinear Maxwell theories with well-defined 
Galilean limits (and correspondingly generalized Yang-Mills equa- 
tions), and a linear modification motivated by the coupling of the 
electromagnetic potential with a certain nonlinear Schrodinger 
equation. In the final section, revisiting an old idea of Lorentz, 
we write Maxwell's equations for a theory in which the electro- 
static force of repulsion between like charges differs fundamentally 
in magnitude from the electrostatic force of attraction between 
unlike charges. We elaborate on Lorentz' description by means 
of electric and magnetic field strengths, whose governing equa- 
tions separate into two fully relativistic Maxwell systems — one 
describing ordinary electromagnetism, and the other describing 
a universally attractive or repulsive long-range force. If such a 
force cannot be ruled out a priori by known physical principles, 
its magnitude should be determined or bounded experimentally. 
Were it to exist, interesting possibilities go beyond Lorentz' early 
conjecture of a relation to (Newtonian) gravity. 

It is a pleasure to dedicate this paper to Gerard Emch, whose skeptical per- 
spective helps motivate those who know him to the pursuit of deeper scientific 
understandings. 
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1. Introduction 



Maxwell's equations are among the most beautiful in physics, unifying the 
forces of electricity and magnetism in a classical field theory that explains 
electromagnetic waves [T]. Some well-known, profoundly- motivated varia- 
tions on Maxwell's equations have included the Born-Infeld theory (a non- 
linear but Lorentz-covariant modification, that introduces an effective upper 
bound to the electric field strength), and the Yang- Mills equations (introduc- 
ing non-Abelian gauge potentials) [21 [3] . These ideas go back many decades, 
and have deeply influenced the development of theoretical physics. Indeed, 
there has been a recent resurgence of interest in non-Abelian Born-Infeld 
Lagrangians [4], which turn out to have important application in string the- 
ory and related subjects [SJ El 0, E] ■ More recently, variations of Maxwell's 
equations have been considered as "test theories," with respect to which ob- 
servations in astrophysics can provide upper bounds to deviations from the 
usual equations or laws of physics [91 [10]. We nevertheless seek to approach 
the idea of modifying Maxwell's equations in new ways with appropriate hu- 
mility. None of the variations considered in this article is ad hoc. Rather, 
each occurred in answer to a specific question in fundamental physics. 

Sections 2 and 3 review two such modifications considered by the sec- 
ond author in recent years. The first of these, proposed in joint work with 
Vladimir Shtelen [TTJ [T2] , is a class of Galilean nonlinear Maxwell theories, 
together with non-Abelian versions that generalize the Yang-Mills equations 
and the non-Abelian Born-Infeld equations. These possibilities arise in an- 
swer to the question of whether and how Maxwell's equations for the four 
fields E , B , D and H can survive unchanged in the Galilean limit of 
c —>■ oo , a feature that is present in neither the usual, linear Maxwell theory, 
nor the usual Born-Infeld theory. The second variation is a simple, linear 
modification, that can be associated with a change over time in some "con- 
stants" of electromagnetism. This possibility comes up in answer to the 
question of how to write gauge-invariant expressions for the electromagnetic 
field strengths E and B , when Maxwell's equations are coupled with a 
natural and very general family of nonlinear Schrodinger time-evolutions in 
quantum mechanics [13] . Such a family of nonlinear Schrodinger equations 
was developed and studied in a series of articles, in joint work with Heinz- 
Dietrich Doebner and Peter Nattermann [13} [151 IB] • I n t ne general nonlinear 
equation, a certain "frictional" term originally proposed by Kostin [T7j is the 
one requiring a change in Maxwell's equations for the field strengths, after 
nonlinear gauge transformations are taken into account. 
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Finally, in Section 4, we take up a different question. In a 1900 article 
that seems to be not very well-known today, Lorentz explored the idea that 
(Newtonian) gravity could be explained if the electrostatic force of repulsion 
between like charges were smaller in absolute magnitude than the electro- 
static force of attraction between unlike charges [IS] . Setting aside Lorentz' 
conjectured relation to gravity, we want to reopen the possibility of a differ- 
ence in magnitude between these forces. There is then a straightforward and 
elegant description of the situation by means of Maxwell's equations, that 
was partially written down in Lorentz' original article. Under the given hy- 
pothesis, one may introduce new electric and magnetic fields whose governing 
equations separate into two fully relativistic Maxwell systems — one describ- 
ing ordinary electromagnetism, and the other describing an overall attractive 
or repulsive long-range force that couples to an "absolute charge." While the 
latter force might conceivably have something to do with gravity, it is more 
plausible to regard it as an extremely small, but theoretically possible, correc- 
tion to ordinary electromagnetism. Such a correction could, in principle, be 
time-dependent, and serve as a further "test theory" for astrophysical mea- 
surements. Whether attractive or repulsive, a modified electromagnetism 
could be important in modeling the early universe. Unless a known physical 
principle rules out such a force a priori, its magnitude should be regarded as 
an experimental question. Well-known nonlinear and non-Abelian general- 
izations of Maxwell's equations, and their unification with weak interactions, 
could then equally well be constructed from the new equations, opening up 
interesting possibilities. 

In the remainder of this section we establish notation, summarize some 
elementary, familiar background material, and make a few relevant remarks. 
All of our discussions pertain to (3 + l)-dimensional space-time. 

Let us write Maxwell's equations in SI units, as follows [1J: 

VxE = -— , V-B = 0, VxH= — +j, V-D = p, (1) 

where E(x, t) is the electric field, D(x, t) the electric displacement, B(x, t) 
the magnetic induction, and H(x, t) the magnetic field; p(x, t) is the charge 
density, and j (x, t) the electric current density. 

The first pair of Eqs. ([I]) imply that we can write E and B in terms of 
potentials ($, A), 

dA 

B = VxA, E = — — -V$. (2) 

dt v ; 
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The choice of $ and A is not unique. For an arbitrary smooth function 
0(x, t), new potentials A' and $' that are obtained from the gauge trans- 
formation 

dO 

A' = A + V0, $' = $-— (3) 

at 

give just the same fields E and B . Thus E and B are said to be gauge 
invariant. The condition V • B = expresses the nonexistence in nature of 
magnetic monopoles, which in this article we do not consider changing. Our 
choice of SI units here is motivated by the desire to avoid incorporating the 
speed of light c into the definitions of any of the fields, as we shall later be 
interested in considering the c — > oo limit. 

From the second pair of Eqs. ([1]), there follows immediately the equation 
of continuity, 

dp 
dt 

expressing conservation of electric charge. 

Thus far, the system is underdetermined. To complete Maxwell's equa- 
tions in the presence of matter, one introduces constitutive equations relating 
the fields D, H to the fields E, B. The usual system of Maxwell equations 
in vacuo is obtained using linear constitutive equations, 



+ Vj 







(4) 



D = e E, H = — B 

fJ'O 



(5) 



where eo^o — 1/c 2 . However, we shall shortly be considering a certain class 
of nonlinear constitutive equations. 

One may take the point of view, given the absence of magnetic monopoles, 
that the only physically detectable fields are E and B. These are defined 
operationally via the observed Lorentz force F on a small "test particle" 
with electric charge q moving with velocity v: 



F = gE + gvxB 



(6) 



The fields H and D can then be regarded as unobservable constructs used 
to describe, by way of the latter two Maxwell equations and the constitutive 
equations, how the observable fields are produced by charges and currents. 

As noted in Ref . [11] , we then actually have a more general class of linear 
constitutive equations, 



(7) 
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which — when combined with Eqs. ([T]) — lead for all values of the real param- 
eter A to the same set of equations for the observable fields E and B. Thus 
the choice A = resulting in Eqs. §5§ is arbitrary. 

Also relevant to the forthcoming discussion is the well-known minimal 
coupling of the electromagnetic potentials ($, A) with Schrodinger's equa- 
tion for the (complex-valued) wave function ip(x.,t) of a single quantum- 
mechanical particle having charge q, 

dib 1 

ih-£- = i-ifiV - qA) 2 4> + q$?p . (8) 

at 2m 

Local U(l) gauge transformations act on ib according to the formula, 

V>'CM) = e^ (x '*V(x,t), (9) 

and it is easily checked that if ib obeys Eq. (jSJ), then ip' obeys a gauge- 
transformed equation of the same form, with new electromagnetic potentials 
given by Eqs. in which 0(x, t) = (%/q) 9(x,t) . 

The gauge-invariant fields E and B, which exert the electric and mag- 
netic forces on the charged quantum particle, are obtained from $ and A 
using Eqs. (T5J), and satisfy Eqs. ([T]). The gauge-invariant probability and 
probability flux densities for the particle are given, respectively, by 

p gi = ^, J gi = hpVip - CVfytb] - —ipipA; (10) 
2im m 

these also obey an equation of continuity. 

Let us remark on the fact that Eqs. ([1]) respect the Lorentz transfor- 
mations of special relativity, while Schrodinger's equation respects Galilean 
transformations. The minimal coupling of Eq. ([8]) is compatible with these 
facts because Eqs. ([T|) also respect Galilean transformations. It is the linear 
constitutive equations that impose Lorentz symmetry on the usual Maxwell 
equations, breaking the Galilean symmetry (see below). 

We reproduce Lorentz transformations here in SI units for completeness. 
Let the subscript || indicate the component of a vector in the direction of 
the velocity v of an inertial frame of reference, let the subscript _L indicate 
the component perpendicular to v, and let v — |v|. Then with 

7 = 7r^' 



5 



the space-time transformation under the Lorentz boost is 

xj|= 7 (X||-Vt), x' ± =x ± , t' = 7 (*-^) ; (12) 
the field transformations are 

Bj|=B| h B; = 7 (B-ivxE) i , Ej|=E| h E'_l = 7(E + vxB) ± , 
H[|=H| h H'_l = 7(H- v x D) ± , (13) 
D;,=D||, D^= 7 (D + ivxH) i; 

II - c z 

and the electric current and charge density transformations are 

j / H=7(j||-p^), j' ± =j±, p' = 7(p-^). (14) 
The corresponding electromagnetic potential transformations in SI units are 

$' = 7 ($_v A), A' = 7 (A-^-$). (15) 

c 2 

Under Lorentz transformation, the following combinations of the fields are 
then invariant: ( 

h = B 2 E 2 , I 2 = B • E ; 

c 2 

h = D 2 - \ H 2 , I 4 = H D ; (16) 
<r 

J 5 = B H - E D , 7 6 = B D + 4e H. 

c 2 

The Born-Infeld Lagrangian as a function of these invariants is 

£ = !-*, R-JLJi + p,-*!*. (17) 

Next let us write the above in covariant notation. Define x^ = (ct,x), 
\i = 0,1,2,3, and x M = g liv x u = (ct, — x), where the metric tensor g^ u = 
diag(l, —1, —1, —1), and where summation over repeated Lorentz indices is 
understood; for example, x^ = c 2 t 2 — x 2 . We further define <9 M = d/dx^ = 
[ (1/c) d/dt, V ]; and we shall use the antisymmetric Levi-Civita tensor e"^", 
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with e 0123 = 1. Then the usual relativistic tensor fields F a/3 and JF Q/3 , con- 
structed from the fields E and B, are 



a/3 



a/3 



r 

with 
Likewise, 








(1/c)^ ( 


l/c)E 2 (l/c)E 3 ~ 




(l/c)^ 





—B 3 B 2 




(l/c)E 2 
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-B 1 
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)E 3 
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B 1 









-B x 


— B 2 ~B 3 




Bi 





a/c)E 3 -(1/ 






B 2 


-{l/c)E 3 


(l/c)£! 




B 3 


a/c)E 2 


-{l/c)E x 


= 0°V 


Efj,l/ 1 J~< 0lj3 


— goi^g^J- • 





(18) 



(19) 
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, C Q/ g — g a fj,gi3uG^ v , (20) 



21) 



and so forth. Maxwell's equations ([T]) then become 

with j 13 = (cp, j). With = ($, — A), we have from the first of Eqs. (1211 . 

j^ = ^d u A v , F ia , = d lt A v -d v A». (22) 



The first two invariants of Eqs. f|T6|) . that enter Eq. (fTTj) for the Born-Infeld 
Lagrangian density, are now written 



h 



2 P v ' 



._ TP t^ u 



(23) 



Note that our strategy, following Refs. [TT] and [12], has been to postpone 
writing constitutive equations for as long as possible. These now relate D, H 
to E, B; or, equivalently, they relate G a @ to . The general form for 
Lorentz-invariant constitutive equations is given by [19] 



D = MB + —NE , 1 1 NB-ME 



(24) 



where M and N are functions of the Lorentz invariants in Eqs. (1161) . or 

B = KD + — QU, E = QD - RH, (25) 
cr 

where Q and R are likewise functions of the invariants. The linear constitu- 
tive equations (jSj) correspond to the choices M — 0, N = 1/fAo, in which case 
eo = 1/ jiQC 2 . Other choices lead to nonlinear relativistic field equations, such 
as Born-Infeld or Euler-Kockel electrodynamics. It is natural in Eqs. (1241) to 
take M and N to be functions of just the first two invariants l\ and I 2 
(which depend only on E and B ), or in the inverted Eqs. (1251) to take R 
and Q to be functions of just the invariants J 3 and 7 4 (which depend only 
on D and H ) . 

The usual approach to writing nonlinear Maxwell theories is to begin 
with the Lagrangian £ , which can be written as a function of the invariants. 
The constitutive equations then follow from the Euler-Lagrange equations. 
But not all Lorentz-covariant theories are Lagrangian, and the approach via 
constitutive equations is more general. 

In covariant form, Eqs. fl2^|) become 

B T BT 

= NF^ V + cMT^ = Mi— — — + M 2 - — — , (26) 

B F^ u BF^ V 

where Mi and M 2 are likewise functions of the Lorentz invariants. 

Now we are ready to discuss the three variations to which this article is 
devoted. 



2. Maxwell equations having Galilean limits 

When Maxwell's equations are written for the four fields E, B, D, and H 
as in Eqs. ([T]), the system is underdetermined. In the SI units we are using, 
these equations are independent of the speed of light c. Furthermore, the 
corresponding equations for the primed fields defined by Eqs. ffl3l) and the 
primed currents defined by Eqs. fTHj) . when written in the primed coordinates 
defined by Eqs. f fl2|) . are unchanged from Eqs. ([T|) — even though the Lorentz 
transformations given in Eqs. f|T2|) - fri4|) are parameterized explicitly by c. 
Since the invariance of Eqs. ([T]) holds for every value of c, we should not be 
at all surprised that in the limit c — » 00, these equations also respect the 
resulting Galilean transformations, 

X ' = x — vt , t' = t , 
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B' = B, E' = E + vxB, 



H' = H-vxD, D=D, (27) 
j'=j-pv, p' = p. 

Now, the value of c does appear in the constitutive equations. Thus the 
choice between Lorentz or Galilei symmetry, or the selection of a particular, 
finite value of c under which the Lorentz symmetry holds, resides entirely in 
the constitutive equations. 

For the linear constitutive equations (JSJ) or (J7J), the speed of light is spec- 
ified by c = (eoHo)' 1 / 2 . In this case, taking a Galilean limit requires that 
some aspect of Maxwell's equations be sacrificed, as discussed in detail by 
Le Bellac and Levy-Leblond [2D] . 

The constitutive equations that select Galilean symmetry, when combined 
with Eqs. ([I]), are 

D = MB, H = NB — ME, (28) 

or, equivalently 

B = RD, E = QT> - RH, (29) 
where M and N, Q and R are arbitrary functions of Galilean invariants, 

h = B 2 , J 2 = B ■ E; 

h = D 2 , h = H D; (30) 

J 5 = B H - E D, J 6 = B • D. 

These constitutive equations and field invariants are respectively the formal 
limits as c — > oo of their Lorentz invariant counterparts. 

As discussed in Ref. [20], however, taking the mathematical step of letting 
c — ► oo is not precisely the same thing as imposing the low velocity condition 
v/c « lona class of physical systems governed by the dynamical equations 
with Lorentz symmetry. For instance, when E and B are held fixed, the limit 
c — ► oo does not allow the "electric limit" of Ref. [20], although v/c « 1 is 
compatible with it. 

Letting M be a constant in Eqs. (1281) requires (since V • B = 0) that 
the charge density p = 0. Hence the answer to the question of a consis- 
tent Galilean electrodynamics, retaining Maxwell's equations, the continu- 
ity equation, and the Lorentz force, is a class of essentially nonlinear the- 
ories, that can arise as the c — > oo limit of a class of essentially nonlinear 
Lorentz-covariant theories. Indeed, Le Bellac and Levy-LeBlond emphasize 
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(always assuming linear constitutive equations) the mutual incompatibility 
of Galilean invariance, the continuity equation with non-zero values, and 
magnetic forces between electric currents. 

With nonlinear constitutive equations, these features are no longer in- 
compatible. Nontrivial choices of M and N in Eqs. (1281) . or R and Q in 
Eqs. (1291) . combined with Maxwell's equations, yield fully consistent Galilean 
versions of electrodynamics. 

For example, Ref. [11] proposes to set 

Q = - R = a + 2a 2 e , (31) 

e \Tj\ 2 

which are homogeneous functions of the field strengths. This can be shown 
to lead to an interesting, albeit non-Lagrangian, theory. 

Writing the Lagrangian for a nonlinear relativistic theory as 

C = C(h,I 2 ), (32) 

a short calculation in Ref. [12] demonstrates from the Euler-Lagrange equa- 
tions that 

N=2 w M =-|- (33) 

Therefore the necessary compatibility condition for the constitutive equations 
to describe such a Lagrangian theory is given by, 

In the Galilean limit, we would take C = C{I\, I2), and argue similarly. 

The usual Born-Infeld theory does not have a nontrivial Galilean limit. 
The Lagrangian £(I\, I2) is given by 



c 2 c 2 



c = 1 ~ n ' n = J^T + v h -¥'l' < 35) 

which leads to the constitutive equations (12^1) with 

M(/ -« = ^fe' N =^R- (36) 

Taking c — >• 00, we have of course I\ — > I\ = B 2 , and I2 — > I2 = B • E . But 
for large c, one has 1Z ~ (c/6)[/i — 7|/6 2 ] 1//2 , whence the limits of M and N 
are both zero. 



10 



It is therefore suggested in Ref. [12] to modify the Born-Infeld Lagrangian, 
replacing 1Z by 




where Ai, A2 are new constants with the dimensions of 1/c 2 . Now, taking 
c — ► 00 , one obtains the Galilean constitutive equations ff28|) . with 

M= X2h N = bXl (38) 

poby/Xji - X 2 I 2 2 /b 2 (JoyJXJi - X 2 I 2 2 /b 2 

Similarly, generalizations of classical (non-Abelian) Yang-Mills theory 
are written by means of Lorentz-covariant, nonlinear constitutive equations. 
Again, with appropriate choices for the dynamics, the new systems can have 
fully consistent Galilean-covariant limits as c — > 00. In analogy with Eqs. 
(I37p - (j38p . one obtains a class of generalizations of non-Abelian Born-Infeld 
theories that are of this type [T2] . 

3. Modification from a nonlinear Schrodinger equation 

Another variation on Maxwell's equations occurs as a result of considering 
the coupling of external electromagnetic fields with nonlinear Schrodinger 
time-evolutions [T31 US]. First we write the class of Schrodinger equations 
under consideration. Refs. [H| and [TJ] provide extensive motivation and 
development, that we omit here; we mainly follow the discussion in Ref. [T3] . 
Letting ^(x, t) be the quantum-mechanical wave function, and define 

p(x,t) = j(x,t) = (l/2z)[^W>-(W#]. (39) 

In this article we shall use the notation p (x, t) to refer to the spatial proba- 
bility density for the quantum-mechanical particle, to distinguish it from the 
net charge density p (x, t) that appears in Sees. 1 and 2. Thus p is here the 
expression we called p gl in the first of Eqs. (|10|) . 

Define the real, homogeneous functionals R\ [ip] , . . . , R5 [ip] , by 

R\ = — i-^j R2 = — R3 = 4j, Ri = - ^ - , R$ = - !*) ■ (40) 

p p p z p z p z 
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The Laplacian in the linear Schrodinger equation (jSJ) can be expanded with 
respect to this basis of functionals, 



(41) 



so that it does not appear explicitly in the equation we shall next write down. 
The general family of nonlinear Schrodinger equations takes the form, 



L j=l 



P 



3=1 



V- (A(x,t)/3) A 2 (x,t) j 
+ Uht,t) H 1 h «i hip + a 2 S, (42) 

P P 

where ip is dip/dt, the coefficients z/,- (j = 1,2), \ij (j = 1,... ,5), and 
atj (j = 1,2) are all continuously differentiable, real- valued functions of t ; 
S(x,t) is arg[^(x, t)] ; U(x,t) is a real-valued scalar function; and A, Ai, 
and A are distinct real-valued, time-dependent vector fields. 

Writing a class of nonlinear Schrodinger equation by adding terms of the 
form (HOI) to the usual, linear Schrodinger equation as in Ref. [I5j, we have 



di/j 



H ip + -KD R 2 [ip]ip + h^, Dj Rj 

3=1 



1>- 



(43) 



and D and the D- 
is obtained from 



where H ip is given by the right-hand side of Eq. ([HI 
have the dimension of diffusion coefficients. Then Eq. 
Eq. (T42|) with the values, 

h l n < Q . 

v\ = , v 2 = i A = A , 

2m 2 2m 

i h i ^ , , h , 

Pi = D 1 , /i 2 = -— + D 2 , p 3 = — + L> 3 , /i 4 = J D 4 , p 5 = — + D 5 
4m 2m 8m 



a i 



A 2 , Ai = 0, A 2 



A 



m 



«2 



0. 



(44) 



The motivation for the form adopted in writing Eq. (|42|) . for the pres- 
ence of the terms with a.\ , a 2 and A\ ^ 0, and the time-dependence of the 
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coefficients, is the possibility of introducing a group of nonlinear gauge trans- 
formations that leave this family of equations invariant (as a class). With 
ip = Rexp [iS], these take the form if) \— ► if)' = R' exp [i<S"], with 

R' = R, S' = AS + ^\nR + 6; (45) 

where 7 and A are continuously different iable, real-valued functions of t , 
A ^ 0, and 6 is a continuously differentiable, real-valued function of x and 
t. Then nonlinear gauge transformations obey the group law, 

(Ai, 71, 0i ) (A 2 , 72, 02 ) = (A1A2, 7i + A l72 , 0! + K x e 2 ) . (46) 

With quantum-mechanical measurements characterized as sequences of po- 
sitional measurements (at distinct times), together with the application of 
external fields between positional measurements [211 [22] . and maintaining 
the standard interpretation of p — \if)\ 2 as a probability density for the 
outcomes of positional measurements, such transformations then leave the 
distribution of outcomes of all measurements invariant. They are also local 
in space-time, and respect a separation condition for multiparticle product 
wave functions [16, 23]. Notice that with 7 = and A = 1 , we recover the 
usual local U(l) gauge group of electromagnetism (acting linearly on if>) as 
a subgroup of the larger group of nonlinear gauge transformations. 
Under the nonlinear transformation in Eq. (14"5]) . we have 

p> = ^'^> = p, 

j' = = Aj + \Vp + pV0; (47) 

so that p (as desired) is gauge-invariant (for nonlinear as well as linear 
gauge transformations). Thus we shall also write p = p gl when we want 
to emphasize this. But j is not gauge-invariant — one must write a new 
gauge-invariant current (see below), to fully generalize Eqs. ([TO]) . Moreover, 
if if} satisfies Eq. (j4"2l then ip' likewise satisfies an equation of this form, 
but with new (primed) coefficients and new external fields. The coefficients 
transformed under (7, A, 9) are given by 



ZA = -r 



A 



^' = -^1 + ^2, 



/^i = -^i + Mi, f4 = 7^"! ~ 7^2 - -Pi + Ap 2 
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/ _ i 7 , / _ 7 2 7 i a 

• 7 1 ( A \ . A . , 

% = Aot! - - a a + - I — 7 - 7 I , a 2 = a 2 - . (48) 

Observe that even if one begins with ol\ = a 2 = and with time-independent 
coefficients z/j and fij , the fact that 7 and A can be time-dependent requires 
that in Eq. f|4*2l) the a, be permitted to take nonzero values, and that all 
the Uj, fij , and otj be, in general, time- dependent. However, the nonlinear 
gauge transformations do not mix the coefficients otj with the coefficients 
Uj and fij . 

The nonlinear term with coefficient a.\ was first proposed as a modifi- 
cation of linear quantum mechanics by Bialynicki-Birula and M icielski [24J , 
and the term with coefficient a 2 was proposed still earlier by Kostin |17j . 
The term with coefficient v 2 7^ was derived by Doebner and Goldin from 
considerations of local current algebra representations [H] , which led to the 
generalized equation containing the nonlinear functionals Rj |15j . 

The external fields of Eq. (142]) transformed under (7, A, 9) are given in 
Ref. Q3 by 

A' = A - ^W, 
A 

A[ = AAi - iA - \M+ (v^i- + ^3-^4) V0, 

Ai = A 2 -^V9, 
U> = AU _ + ^_ a2 j + ^[W] 2 + 

(//4-/X3 J) V 2 £ + |V-A - A • V0. (49) 



Observe that Eqs. fj48|) - fj49|) imply both the nonzero A\ and nontrivial A 2 
values are required in Eq. (|42p . Even if one begins with A\ = and 
^4 2 = —2.4., as in the linear Schrodinger equation [c/. Eqs. dSJ)], the 
nonlinear gauge transformations compel one to introduce more general values 
for these fields. Nonlinear Schrodinger equations with arbitrary values of A 2 
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were proposed by Haag and Bannier [25] , while the interaction with a general 
external vector field A\ was considered in Ref. |2bj . 

Next let us write the equations of motion described by this class of non- 
linear Schrodinger equations entirely in terms of gauge-invariant quantities, 
as in Ref. [13] — where "gauge invariance" is interpreted with respect to the 
group of nonlinear gauge transformations. To start, a gauge-invariant cur- 
rent density J gl may be written 



2i/ij - 2z/ 2 Vp - 2pA, (50) 



which evidently reduces to Eq. fllUp for the linear Schrodinger equation [when 
v\ = —h/2m, v 2 = 0, A = (g/2m)A]. We thus have (again) a continuity 
equation for the probability density, dp gl /dt = — V ■ J gl . Refs. [15] and 
[T6] provide a set of gauge-invariant parameters necessary for the desired 
description, 

Ti = v 2 -~Hi, t 2 = VxlL-i - villi , r 3 = — , r 4 = /i 4 - / ui — , 
2 v\ V\ 

, 2 ^3 
T 5 = Vi[L h - z/ 2 /i 4 + v 2 — ) 

fix = v x ai - v 2 "2 + v 2 — - i>2 , (3 2 = a 2 - — ■ (51) 

When T\ ^ 0, T4 ^ 0, or f3 2 7^ 0, time-reversal invariance is violated. When 
T3 7^ —1 or T4 ^ 0, Galilean invariance is violated. 

Gauge-invariant external magnetic and electric fields are now given by 
the formulas [T3] 

B = V xA = ^-B, 

2m 

£ = -™- d 4-^ A = £ E > (52) 

where 

U = - u t U - r 3 A 2 - (r 4 - 2nr 3 ) V • A + A ■ A 2 - v* V • A 2 ■ (53) 

Thus U is to be identified with (g/2m)$ (which may be directly checked for 
the linear Schrodinger equation); but the main point here is that the formula 
for E in terms of $ and A has been modified from Eq. (J2]) to include an 
extra term, originating with Kostin's nonlinearity: 

E = -V$ - — - (3 2 A. (54) 
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The extra term is necessary — if we leave it out of Eq. (I52p . £ fails to be 
gauge-invariant. We also have new gauge-invariant external vector fields, 



( 2z/ 2 /i 3 

V "I 



/ii — n± 1 A — V2A2 , A 2 



2p 3 



A 2 - A . (55) 



Continuing to follow Ref . [13] , we are now in a position to write equations 
of motion obtained from Eq. (jl2"j) . in a form that is manifestly gauge- invariant 
for the group of nonlinear gauge transformations. We use the hydrodynamical 
variables variables p gl and J gJ /p gl , where the latter has the interpretation 
of a gauge-invariant velocity field. We then have the interpretation of the 
(gauge-invariant) magnetic induction as a vorticity field, 



V x 



J gi 

~P^ 



-2B 



m 



B 



(56) 



The dynamical equations are the equation of continuity that we wrote above, 

<9p gi 



dt 



-V ■ J gl 



(57) 



together with the equation for the changing velocity field, 



d_ 

dt 



jgi 



V 



2n V 



jgi 

p^ 



2r 2 



V 2 p gi 



r 3 



jgi 

p^ 



+ V 



+ V 



2n [1 + r 3 ] - r 4 



V- (Afp^ 



P 



Jgi 

P^ 



2t 3 A? 



P 



gi 



jgi 

p^ 



2/3i lnp 



gi 



p gl / m 



(58) 



In Eq. (1581) we see that /5 2 , taken to be positive, has a natural interpretation 
as a gauge-invariant coefficient of friction — it governs the magnitude of the 
term in d t (J gl /p gl ) that is proportional to J gl /p gl . 

In the framework of the nonlinear quantum mechanics discussed here, the 
(gauge-invariant) expected values for the position, velocity, and acceleration 
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of the quantum particle (all of which are functions of t) are given respectively 
by the following expressions: 



< x > 



x p gl (x) (ix . 



< V > 



d < x > 
dt 



iS 1 



Jgi 

r ] dx 

P gl 



J gi (x)dx. 



(59) 



< a > 

2 

+ 



(9 < v > 



jgi 



x 

m 



d_ 

dt 



jgi 

P^ 



c?x . 



In Eqs. (I58l) -( l59l . we see that the laws of force describing the interaction of 
the charged particle with the E and B fields are unchanged from those in 
linear quantum mechanics. 

Now Eq. fl5^|) gives us the variation on Maxwell's equations that is the 
focus of this section. The usual equations for E and B are replaced by 



V x E 



OB 

~dt 



- B , V • B = . 



(60) 



Note that the second of these equations is still consistent with the first. Let 
us take 02 to be a constant, independent of t . If E (x, t) , B (x, t) satisfy 
the original Maxwell equations (with (3 2 = 0), then fields satisfying Eqs. (!60l) 
are given by 

E = E e~ f32t , B = B e~ l32t . (61) 

But nonlinear quantum mechanics alone does not specify the remaining two 
Maxwell equations. One possibility is to maintain the constitutive equations 
([5]), with fixed coefficients Eo and /io . Then taking E , B , D = £oEo , 
H = (l//i )B , p and jo to satisfy the usual, linear Maxwell equations 
([T]), we must have 



D = D e 



-Pit 



H = H e 



and 



V ■ D = p , with p 
<9D 



Poe 



-fat 



(62) 



VxH = — + /3 2 D+j, with j = jo e 



-Pit 



(63) 
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That is, with (3 2 > , the magnitudes of all the electric charges and currents 
are decaying exponentially with time. Of course, the equation of continuity 
for p and j no longer holds, and net charge is no longer conserved. Instead, 
we have the equation 

^ = -V-j-/3 2P . (64) 

The nonzero value for /?2 , interpreted as a coefficient of friction experienced 
by a charged particle, has introduced a preferred universal reference frame. 
The corresponding Maxwell theory is no longer covariant. 

Another alternative is to join the standard Maxwell equations for D 
and H with Eqs. (1601) by supposing that the coefficients in the constitu- 
tive equations — the permittivity and permeability of free space — are time- 
dependent. With 

e{t) = e e + ^\ pit) = floe-**, (65) 
we may combine Eqs. (I6"U|) and fIBll with the constitutive equations, 

D = e(t)E, H=^-B. (66) 
p{t) 

Then D = D , H = H , p = p , and j = j . We have e(t)p(t) = 1/c 2 
for all t . The equation of continuity holds for p and j , and the net charge 
is conserved. 

Since current systems of units define the vacuum permittivity to have 
a fixed numerical value, the modification embodied in Eqs. (166]) requires 
some reexamination of the way in which we define our units of measurement 
for electromagnetism. We shall return briefly to the idea of time-dependent 
permittivity and permeability toward the end of the next section. 



4. Force differences between like and unlike charges 

For our final variation on Maxwell's equations, we pursue the suggestion 
by Lorentz that the electrostatic force of repulsion between like charges is 
slightly different in magnitude from the electrostatic force of attraction be- 
tween unlike charges. While Lorentz proposed this idea in 1900 as a way to 
arrive at a universally attractive (Newtonian) gravitational force, we want 
to regard it here as just a modification of classical electrodynamics. Let us 
write a description of the situation in terms of field strengths. 
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The electrostatic forces among pairs of positively-charged point particles 
having charges Q + > and q + > and/or negatively- charged point parti- 
cles Q_ < and g_ < , is described by 



1 Q+Q+ 



1 Q_g_ 



, , or (67) 

where r is the unit vector at the location of each particle, pointing away from 
the other, and r is the distance between the pair. The subscripts I and u 
stand for "like" and "unlike" respectively, and Eg ^ e u . If \Q+\ = \Q-\ and 
|<7+| = \Q-\i then the forces F u in the second equation are equal. 

Introduce the electric field E + exerting force on positive charges, and the 
electric field E_ exerting force on negative charges; also the displacement 
fields D + and D_ produced (respectively) by positive and negative charges. 
In ordinary electromagnetism, E = E + = E_ , while D = D + + D_ . 
Evidently we must keep track separately of the density p + > of positive 
charge and the density p_ < of negative charge. From Eqs. (1671) . the 
force F + experienced by the positively charged particle with charge q + in 
the presence of a composite having charges Q + and Q- at a distance r is 
given by 

+ 1 Q- 

47T££ r 2 4:7re u r 2 



g+E H 



(68) 



and similarly for a negatively charge particle, 



g_E_ 



Q- 



4:7re u r 2 



+ 



1 Q- 
Attei r 2 



(69) 



There is a net electrostatic force between neutral composites that is attractive 
if e u < Ei , and repulsive if e u > Ei . The situation thus far is described by 
the Maxwell equation 



and the constitutive equation 







P+ 


D 




. P- _ 









l/e« 




D + 


E 




l/e« 









(70) 



(71) 
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Similarly, introduce B + and B_ as the magnetic inductions that exert 
velocity-dependent forces on positive and negative moving charges q + and g_ 
(respectively); so that the total forces are given (respectively) by 



Q 



_v x B. 



g_E_ + g_v x B_ 



(72) 



Evidently, we must also keep track separately of the current of positive charge 
j + and the current of negative charge j_. We introduce H + and H_ as the 
magnetic fields produced (respectively) by electric currents j + and j_, and 
(respectively) by changing displacement fields D + and D_ . Then we obtain, 
consistent with Lorentz covariance and the absence of magnetic monopoles, 
the additional Maxwell equations, 



V x 



H 
H 



d_ 

dt 



D 
D 



+ 



V x 



with the constitutive equation 





d 


B + " 


, v- 


B + 


E 




B 


B 



o, 







Hi 


Hu 




H + 


B 




_ Hu 


Pi _ 




H_ 



where 



1 

£lHl — £ uHu — ~t; 



(73) 



(74) 



(75) 



From Eq. (1701) and the first of Eqs. (173]) . we have separate continuity equa- 
tions for p+ , j + and p_ , j , 



d_ 
dt 



p+ 


+ V- 


j + 


= 


. p- . 




. j- . 





(76) 



Thus far we have written in Eqs. (I70|) and (1731 a doubled set of Maxwell 
equations, indexed by (+, — ) , that are coupled by the matrix constitutive 
equations (j7T]) and (I74|) . To recover the theory as a perturbation of ordinary 
classical electrodynamics, define p = p + + p_ and D = D + + D_ ; so that 
p is the net charge density, and V ■ D = p. Defining E = |(E+ + E_), we 
recover the constitutive equation E = (l/e Q )D by setting 



1 

so 



2 \e t 



1 



(77) 
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Similarly, defining j = j + + j and H = H + + H_ , we have the usual 
Maxwell equation V x H = dD/dt + j. Letting B = |(B+ + B_), we 
also have the Maxwell equations V x E = —d~B/dt and V • B = 0, and we 
recover the constitutive equation B = yU H by setting 

to = ~ [to + to) • (78) 

Then it follows from Eqs. (ED), fTTTTi . and jTSJ that e to = 1/c 2 . In short, 
we still have the four fields E , B , D , and H , and they still satisfy the 
usual Maxwell equations incorporating the net charge density and net electric 
current density, with coefficients e and /i = 1/ e c 2 . 

But we also have an additional set of fields, density, and current, 

D = D + — D , E = i(E+-E_), 

H = H + - H , B = l - (B + - B ) , (79) 

P = P+- P-, j = J+ ~ J • 

These also obey Maxwell's equations; but with new constants in their con- 
stitutive equations, e and Jl , that are given by 

1 1/1 1\ _ 1, 

~ = « - - — ), V =7i (to -to)- (80) 
e i \£i e u j z 

In ordinary electromagnetism, e is infinite and Jl is zero — so that, although 
p , j , D and H are defined and nontrivial, E = and B = . 

If it is not actually infinite, the magnitude of e is presumably very large 
compared with that of e ; while if it is not zero, the magnitude of Jl is 
small compared with that of /io . One verifies straightforwardly that 7J1 = 
1/c 2 ; but the new constants may be of either sign — both positive (for a net 
attractive force between neutral composites), or both negative (for a net 
repulsive force). We thus have obtained a pair of fully decoupled Maxwell 
systems, consistent with the equations set down in Lorentz' original article. 

An idealized composite point particle, having positive charge q + > and 
negative charge q_ < , may be equivalently described as having net charge 
q = q+ + q~ and absolute charge q = q + — g_ > . Such a particle, moving 
with velocity v , experiences according to Eqs. (1721 the total force 

F = F+ + F_ = q (E + v x B) + q(E + v x B) . (81) 
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Thus the new fields E , B , D , and H couple to the absolute charge (which 
is always positive), and the absolute current; while the usual fields E , B , D , 
and H still exist in this framework and couple to the net charge and the net 
current. One should no longer automatically take positive charge flowing to 
the right to be indistinguishable mathematically or physically from the same 
amount of negative charge flowing to the left. The net currents j are the 
same in these two cases, but the absolute currents j are equal and opposite. 

Let us close this section with some comments and speculations, many of 
them rather obvious, about such a "doubled electromagnetism" theory. 

In his original article [T5], Lorentz took the new force to be attractive. 
He sought to identify the absolute charge with mass (and consequently the 
absolute current with momentum), and to calculate whether the precession 
of the perihelion of Mercury's orbit could then be understood as due to the 
(very small) magnetic force that would originate from the absolute current. 
He concluded that the resulting force would be too weak to explain the 
astronomical observations, and of course this line of thinking was superseded 
by the success of Einstein's general relativity. 

However, we want to entertain the idea of a modified electrodynamics that 
does not identify the extra fields with gravity, and that treats absolute charge 
not as mass but as an additional property of matter. Of course, this does 
not preclude the possibility that existing measurements of gravitational forces 
have erroneously incorporated a small extra electrostatic force (attractive or 
repulsive). It would seem to be an especially interesting conjecture that the 
new force is repulsive. In any case, we appear to have an additional parameter 
with which to fit cosmological models, and an additional "test theory" for 
study through observations in astrophysics. 

Since we have a new set of fields obeying Maxwell's equations, we would 
also need to have a new type of electromagnetic wave (coupling weakly with 
the absolute charge), a new type of photon, and a new quantum electrody- 
namics (see below). 

Now the absolute electric charge of a system must be at least equal to the 
net electric charge, and at least equal to the sum of the absolute charges of the 
system's components. However, it could in principle be greater. The absolute 
charge of a nucleon, for example, might be the sum of the absolute values of 
the charges of its constituent quarks; but one could also conjecture additional, 
unobserved positive and negative charges in equal measure, contributing to a 
larger value of the overall absolute charge. While net charge is quantized in 
fixed units, it is plausible but not necessary that absolute charge be similarly 
quantized. Thus, it does not appear to be inconsistent to take the absolute 
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charge to be proportional to the mass, as Lorentz implicitly did. 

However, modern particle physics offers no fundamental theoretical reason 
to make such an assumption. If we make reference only to constituent quarks 
and leptons, the absolute charge of a proton (comprised of two up quarks and 
one down quark) is 5/3; that of a neutron (comprised of one up quark and two 
down quarks) is 4/3; and that of an electron (taken to be fundamental) is 1. 
Then the absolute charge of a proton together with an electron is 8/3, double 
that of the neutron, while the respective masses are very close to equal. Un- 
der these assumptions, the absolute charge per gram of electrically neutral 
matter comprised of heavier elements is macroscopically different from that 
of matter comprised of lighter elements, and their accelerations under the 
Earth's absolute electric field would be different in magnitude. Since such 
differences are not observed, we should take the empirically-determined mag- 
nitude of any new inverse- square- law force of "absolute electromagnetism" 
to be small compared with Newtonian gravity — making it extremely small in 
comparison with ordinary electromagnetism. 

Still, one may conjecture that physical "constants" are not actually con- 
stant, but change as the universe ages. We discussed in Sec. 3 the possibility 
of the vacuum permittivity and permeability changing exponentially with t , 
producing a modification in Maxwell's equations and a universal frictional 
force that breaks covariance. One may instead take these to be fixed at Eo 
and no , but entertain the possibility that e u and ee are changing, and that 
they have not always been as close as they are today. For example, we could 
have 

e(t) =e e + P\ fi(t) = fi , (82) 

and modify the Maxwell equations for the perturbing fields E , B , D , and 
H as in Sec. 3. 

Furthermore, should there be regions of space-time containing plasmas 
of electrons and positrons, these might contribute proportionally more to 
absolute charge than to gravitational mass. Such speculations leave open 
some possibilities for observable effects in astrophysics, even for a small force. 

While net charge and absolute charge are both conserved when there is 
no particle creation or annihilation, it seems clear that absolute charge is 
not conserved by fundamental particle processes. [Of course, at the time of 
Lorentz' paper, mass, positive charge, and negative charge would all have 
been taken as separately conserved.] As long as we stay with Lorentz' idea 
that absolute charge is proportional to mass, then (as mass is transformed 
into energy during fundamental particle processes) we have a conservation 
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law. But if we take absolute charge to be an independent quantity with 
which the new electromagnetic fields couple, it becomes an unlikely, specu- 
lative possibility that absolute charge transforms into something previously 
unknown during annihilation processes, so as to maintain a conservation law. 
More likely, one should take the new U(l) gauge symmetry to be broken out- 
side the classical domain that is governed by Maxwell's equations, requiring 
a different quantum electrodynamics for absolute electromagnetism. 

We have not discussed the question of absolute charge from virtual particle- 
antiparticle pairs, or vacuum polarization. It appears that, unlike the situa- 
tion for net charge, it should be possible to have a coherent superposition of 
quantum states having different absolute charges. For example, the neutral 
pion, written as a linear combination of up and down quark- ant iquark pairs, 
(uu — dd)/y/2, would combine states of absolute charges 4/3 and 2/3 (in 
units of the electron's charge). 

If Lorentz' conjecture breaks no known physical principle, then the ques- 
tion of a discrepancy in magnitude between the electrostatic forces between 
like and unlike charges is purely an empirical one, and the best we can do in 
principle is to establish an experimental upper bound to this discrepancy (or, 
equivalently, to the ratio e /e). For example, it is clear that his conjecture 
continues to respect the charge conjugation invariance of electromagnetism. 
However, it should be noted that current physics tends to assume the equal- 
ity in magnitude between these forces. Thus the present, official definition 
of the coulomb is effectively as a unit of net charge, while the value of e 
is not measured but fixed by definition. Just as we have come to distin- 
guish (theoretically) gravitational mass from inertial mass, and to regard 
their proportionality as a question to be determined by experiment, so may 
it be necessary to distinguish (theoretically) eg from e u , and to regard their 
closeness as having a value to be bounded by experiment. 

Possible further development of a "doubled electromagnetism" theory 
could entail its involvement in electroweak unification or in the standard 
model, its nonlinear modifications (as discussed for Maxwell's equations in 
Sec. 2 of this article), its coupling with linear and nonlinear Schrodinger 
quantum mechanics (as discussed for Maxwell fields in Sec. 3 of this article), 
its non-Abelian generalizations, and its quantum electrodynamics. 
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